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INTHXDCTXCtr 


Progress by investigators at the Georgia Institute of Technology in the 
development of techniques for passive microwave retrieval of water vapor 
and precipitation parameters using millimeter- and sub-millimeter 
wavelength channels is reviewed. Channels of particular interest are in the 
tropospheric transmission windows at 90, 150, 220, and 340 GHz and centered 
around the water vapor lines at 183 and 325 GHz. Collectively, these 
channels have potential application in high-resolution mapping (e.g. , from 
geosynchronous orbit) , remote sensing of cloud and precipitation 
parameters, and retrieval of water vapor profiles. Both theoretical and 
experimental results to date are discussed. 

The theoretical effort has consisted of radiative transfer modeling 
using the planar-stratified numerical model described in Gasiewski and 
Staelin (1990) . The model incorporates scattering and absorption from 
liquid and frozen spherical Mie-scattering hydrcroeteors. Extremely High 
Frequency (EHF) and Submil 1 imeter-Wave (SMMW) brightness temperatures have 
been calculated for a variety of hydrcmeteor scenarios in order to 
determine the sensitivity of such spectra to hydrcmeteor parameters and 
water vapor. 

The experimented effort has been associated with the design of the NASA 
Goddard Millimeter-Wave Imaging Radiometer (MIR) . This instrument, when 
operated on the NASA ER-2 high-altitude platform, will provide radiometric 
data necessary for evaluation of radiative transfer models and development 
of retrieval techniques. Georgia Tech involvement has included overall MIR 
performance simulation, development of receiver specifications , instrument 
layout, data acquisition system design, and radiometric calibration load 
characterization and design. 
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SOMARJf OF ACTIVITIES 


Three categories of research by Georgia Tech investigators supported 
under this grant can be defined: 1) millimeter- and submil limeter wave 
radiative transfer modelling, 2) performance analysis and design of the 
Millimeter-wave Imaging Radiometer (MIR), and 3) n u m e r ical electromagnetic 
and thermal analysis of microwave blackbody calibration loads. Progress 
within each of these categories will be discussed. 

1. Radiative transfer Modelling 

An analysis of the potential uses of specific Extremely High Frequency 
(EHF: 30-300 GHz) and Sutmillimeter-Wave (SMMW: 300+ GHz) channels for 
passive meteorological remote sensing has been performed using an iterative 
numerical radiative transfer model (Gasiewski and Staelin, 1990) . The 
initial results of this sensitivity analysis are described by Gasiewski 
(1990, see Appendix A) . A more thorough elaboration will be available in 
forthcoming paper to be submitted to the IEEE Transactions on Geoscience 
and Remote Sensing. Reported here are the essential findings: 

A) Clear-air water vapor profiling can be performed, in principle, 
equally well using SMMW channels at 325+\-l,3,7 and 220 GHz as has been 
demonstrated using channels at 183+\-l,3,7, and 166 GHz. The enhanced 
spatial resolution available at the higher frequencies using diffraction 
limited apertures of fixed size make these channels desirable. However, 
water profiling using the 325 and 220 GHz channels alone is expected to be 
more adversely affected by clouds. 

B) Coincident observations at using both lew and high frequency channel 
sets (i.e., 166, 183+\-l,3,7, and 220, 325+\-l,3,7 GHz) are expected to 
provide at least one additional observable degree of freedom useful for 
estimating cloud water content, and possibly, cloud tep altitude. This 
degree of freedom can be seen in the difference between computed 183 and 
325 GHz spectra, and is expected to contain information not available using 
only one of these two channel sets. The additional degree of freedom will 
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be related to cloud particle size and density, but is not expected to be 
roonotonic in these parameters. 

C) Non-precipitating clouds such as high-altitude cumulus and cirrus 
are detectable using a 340 GHz or 325+\- 9 GHz channel. At tropical 
latitudes, the ability to discriminate between water vapor fluctuations and 
non-precipitating clo uds using either of these SMMW channels is somewhat 
degraded relative to that at middle and polar latitu des , although these 
channels are still expected to be useful. E sti m ation of cloud parameters 
such as cirrus ice water content appears feasible, particularly within 
middle and polar latitudes. However, further work cn the retrieval of cloud 
water using EHF and SMMW channels is needed. 

D) The SMMW channels are somewhat more sensitive to clouds than EHF 
channels, although not as sensitive as IR channels. Overall, a 340-GHz 
channel will be less sensitive to storm structure and rain rate as one at 
90-GHz, but more sensitive to thin ice canopies (<0.1g/m3 density) , and 
should be particularly useful for detecting and mapping high-latitude snow 
and ice clouds. However, the sensitivity differences between a 90- and a 
340-GHz channel are net expected to be nearly as significant as between a 
90-GHz and an IR channel. In addition, the increased spatial resolution 
available at 340 GHz relative to 90 GHz (assuming diffraction limited 
apertures of equal size) is expected to benefit storm cell mapping. 

Work is also in progress on the study of the observable modes in 
realistic EHF and SMMW brightness imagery of clcuds and precipitation. 
Hydrameteor profiles are derived frem volume reflectivity data of a 
convective storm observed by the CP-2 radar on July 11, 1986 during 00HMEX. 
The hydremeteor profile data is used as input to the planar-stratified 
numerical radiative transfer model, which subsequently produces high- 
resolution (approximately 1 km spot size) brightness imagery. 

A Karhunen-Ioueve (or principal oenponents) decomposition (Gasiewski 
and Staelin, 1989) is being implemented to reduce the complexity of the 
brightness spectra. This operation allows the sensitivity of the spectra to 
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meteorological parameters to be more easily studied. Using EHF and SMMW 
spectra, observable inodes related to the amount of water and the ice size 
distribution are expected, but have not yet been ascertained. Moreover, 
these modes would be dependent on the spatial resolution of the imagery. To 
simulate satellite-based observations, the high-resolution brightness 
images are to be convolved with appropriate antenna gain patterns, and 
white observation noise will be added. The observable inodes in these 
blurred, noisy images will be representative of what ’can be obtained from 
space-based passive imagers. 

2. Performance A nalysis and Design of the MIR 

The most important activity to date has been the design of the MTR. 
Aircraft data from this instrument will be indispensable for clear-air EHF 
and SMMW radiative transfer and water vapor and precipitation sounding 
studies. Georgia Tech involvement has includ e d overall MIR performance 
simulation, development of receiver specifications, instrument layout, data 
acquisition system design, and radiometric calibration load 
characterization and design. 

Our clear- air radiative transfer calculations suggest that the channel 
specifications in Table 1 should be used to achieve reasonably uniform 
sampling of the vertical water vapor profile. It is noted that the channel 
set 90, 166, arri 183 +\~l/ 3,7 GHz exhibits a series of wei^iting functions 
with more uniformly-spaced peak-altitudes than the channel set 90, 150, and 
183+\-l ,3,7 GHz. Thus, it is expected that the former set (using 166-GHz) 
will be somewhat more useful for water vapor profile retrieval then the 
latter set (using 150-GHz) . In addition, the 166-GHz channel lies in a band 
currently allocated for passive Earth remote sensing. As currently 
envisioned, the MIR will have a channel at 150 GHz (rather than 166 GHz) to 
coincide with a proposed Advanced Microwave Sounding Unit (AMSU) channel. 
This compromise will cause a slight reduction in scientific merit, 
particularly in the ability to experimentally assess any advantages in 
using either 166 or 150 GHz. However, it will allow coincident aircraft and 
AMSU satellite data to be immediately compared. 
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In order to simulate the in-flight performance on board the ER-2 
aircraft, a sp readsheet program was written which allows various instrument 
parameters (e.g. , antenna bea width, receiver noise figure, etc.) to be 
adjusted. The overall effects of these adjustments on many other dependent 
parameters (e.g., scan rate, instrument resolution, data rate, etc.) is 
immediately displayed. This program has been useful in optimizing the 
design of the MIR subject to the constraints of the ER-2 platform, the 
receiver noise figures, the channel frequencies, and. other parameters. As 
currently envisioned, the MIR will have 3.5-degree antenna beawidths for 
each channel. A spreadsheet printout listing this and other dependent and 
independent parameters is shown in Fig. 1. 

The MIR will consist of two units: 1) the scanhead, and 2) the data 
acquisition/control and pcwer system (Fig. 2) . The scanhead will house the 
receivers, calibration loads, scanning mirror and motor, and load 
temperature sensors and regulators, IF and video electronics, and A/D and 
D/A converters. In order to minimize scanhead weight, the MIR will be a 
total power radiometer design. The data acquisition system will consist of 
a computer with floppy drive and cartridge tape mass storage unit. In order 
to retain an acceptable dynamic range and temperature resolution (400 K and 
0.1 K, respectively) in the presence of inevitable receiver gain and offset 
fluctuations, an offset compensation circuit is being designed by Georgia 
Tech. Opto isolation of the analog module will be employed to decrease 
ground-loop noise susceptibility. The power module is to be designed yet, 
but will house all circuit breakers, indicators, DC-to-DC converters, and 
will function as the aircraft interface. 

The implementation of a 325/340 GHz receiver on the MIR has been of 
high priority. We are currently in favor of a double-balanced Schottky 
diode mixer design, with a subharmonical ly-pumped local oscillator (LO) 
Gunn diode source. The ID frequency will be 162.5 GHz; the source will be 
replaceable with a 170 GHz Gunn oscillator for observations at 340 GHz. 
Waveguide, as opposed to quasi-cptical diplexing seems more practical, 
although we are reviewing a single-ended quasi-cptical Schottky design. 
Pending acquisition of funds, we plan to request bids for the 325/340 GHz 
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receiver during March, 1991. The receiver should be ava ila ble for testing 
and integration before 1992. 

3. Numerical Analysis of Microwave Blackb odv Calibration loads 

To accurately calibrate the MIR, the calibration load total 
reflectivity needs to be less than 1%, and known to better than 0.1%. 
Manufacturer's specifications typically provide only the specular component 
of the reflectivity, which is thou^it to be substantially less than the 
total reflectivity. In order to refine our estimates of the total 
reflectivity, we are investigating the scattering of an incident plane wave 
from a calibration load modelled as a lossy periodic wedge grating. 

One approach we are using is close to that described by Maharam and 
Gaylord (1982) , in which the grating is partitioned into a number of thin 
planar slab gratings (see Appendix B) . Within each slab the electromagnetic 
field can be decomposed into Floquet harmonics. A set of coupled wave 
equations is subsequently developed to relate the harmonics in adjacent 
slabs. A com p uter program based on this technique has been implemented, 
although a more economical matrix inversion method is yet needed for 
practical applications. Another approach that we are investigating is based 
on the periodic Green's function (Gasiewski, et ad, 1991) . 

We have also analyzed the steady state thermal distribution in periodic 
wedge and pyramid calibration load structures using thermally conductive 
plate and rod (respectively) models. Initial results suggest that the 
temperature profiles within either wedge- or pyramid-type loads will be 
quite similar (for comparable hei^it-to-pitch ratios) , and hence either of 
these structures are equally acceptable for use as calibration loads in a 
dissipative (i.e. , non-equilibrium) thermal environment. 
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GCNCHEICNS AND HANS FOR FUTURE WORK 


EHF and SMMW radiative transfer calculations have been carried out for 
hydrctneteor— laden tropical , mid-latitude, and polar atmospheres. The 
results suggest that SMMW channels at 325 and 340 GHz offer significant 
potential for passive meteorological remote sensing, and experimental 
confirmation using aircraft observations at these frequencies is urged. 

The main findings have been enumerated in the previous section, and will be 
discussed in detail in a forthcoming publication. 

The Karhunen-Ioueve rank reduction analysis using simulated EHF and 
SMMW brightness imagery provides a quantitative method of determining the 
meteorological information content of passive microwave spectra observed 
using selected channel sets. Future plans include the use of more realistic 
hydrometeor profiles derived from a numerical clcud evolution model (the 
NASA Goddard Cumulus Ensemble - GCE) , as described by Tao et al (1987) . We 
also plan to investigate the effects of an antenna pattern deconvolution 
operator to enhance the resolution of satellite— based imagers. Under 
separate funding, the planar-stratified radiative transfer model is 
currently being upgraded to describe aspherical hydncmeteor distributions. 
This capability will be useful in assessing the sensing potential of 
polarized EHF and SMMW channels for cirrus parameter retrieval. 

To date, we have provided assistance to NASA/GSFC in the mechanical and 
electrical design of the MIR. Georgia Tech's main responsibilities are the 
design of the data acquisition system, implementation of MIR software, and 
development of the 325/340 GHz receiver. Ihe construction of most 
components of the MIR are the responsibility of NASA/GSPC. The earliest 
test flight will be during fall 1991, using only the channels up to 220. 
GHz. The 325/340 GHz channel is expected to be added in 1992. Particular 
attention is being given to the instrument calibration, load temperature 
measurement hardware, video circuitry, and overall system interference 
immunity. 
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Hie radiometric calibration load reflectivity analysis will be 
continued beyond this grant, although this study will be of tertiary 
importance relative to the MIR construction and the theoretical radiative 
transfer analysis. Ultimately, we plan to perform numerical calculations of 
the bistatic scattering coefficient of a lossy periodic grating based on 
coupled-wave and periodic Green's function approaches, and to verify these 
calculations against laboratory measurements using similar wedge-shaped 
calibration loads. The resulting numerical model will be useful for the 
absolute calibration of the MIR as well as in future calibration load 
designs, particularly those wideband loads. Current calibration load 
designs are highly empirical , and the manufacturers' load specifications 
are generally insufficient for accurate calibration. 
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TABLES AND FIGURES 


Table 1. MIR Channel Specifications 


Ch # 

ID Freq 
(GHz) 

XF1 Freq 
(MHz) 

BW 

(MHz) 

1 

89.000 

1000.0 

1000.0 

2 

150.000 

1000.0 

1000.0 * 

3 

183.310 

1000.0 

1000.0 

4 

183.310 

3000.0 

2000.0 

5 

183.310 

7000.0 

2000.0 

6 

220.000 

2000.0 

2000.0. 

7 

325.150 

1000.0 

1000.0 

8 

325.150 

3000.0 

2000.0 

9 

325.150 

7000.0 

2000.0 

10 

340.000 

2000.0 

2000.0 


* Although a similar channel at 166.000 GHz 
would provide a more uniformly distributed set of 
weighting function peaks, 150.000 GHz was chosen 
to coincide with the channel specifications on 
current operational satellites. 
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Millimeter-wave Imaging Radiometer (MIR) 


MIR Data Acquisition System 

GOES Tune 



Figure 2. MIR block diagram, illustrating components in the 
scanhead and data acquisition system. 


13 



